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Abstract 

Here we apply nanomechanical resonators to the study of oscillatory fluid dynamics. A high- 
resonance- frequency nanomechanical resonator generates a rapidly oscillating flow in a surrounding 
gaseous environment; the nature of the flow is studied through the flow-resonator interaction. Over 
the broad frequency and pressure range explored, we observe signs of a transition from Newtonian 
to non-Newtonian flow at lot ~ 1, where r is a properly defined fluid relaxation time. The obtained 
experimental data appears to be in close quantitative agreement with a theory that predicts purely 
elastic fluid response as lot — > oo. 
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The Navier-Stokes equations based upon the Newtonian approximation have been re- 
markably successful over the centuries in formulating solutions for relevant flow problems 
both in bulk and near solid walls [1]. The Newtonian approximation breaks down, how- 
ever, when the particulate nature of the fluid becomes significant to the flow. The Knudsen 
number, Kn = X/L, is one parameter, which is commonly used to settle whether the New- 
tonian approximation can be applied to a medium or not. Here, one compares the mean free 
path A in the medium to an ill-defined characteristic length L. A second defining parame- 
ter, especially for oscillatory flow, is the Weissenberg number, Wi = t/T, which compares 
the characteristic time scale T of the flow with the relaxation time r in the medium. As 
t/T = ujt is varied — for instance, by varying the flow frequency ui or the relaxation time 
t — the nature of the flow changes drastically. 

Recent developments in nanometer scale engineering have created a vibrant subfield of 
fluid dynamics called nanofluidics [2]. Most nanofluidics work is concerned with flow in 
nanoscale channels and remains strictly in the Newtonian regime. In contrast, emerging 
nanometer scale mechanical resonators j^l, 4], with frequencies already extended into the 
microwaves 5|, IfjJ, offer an uncharted parameter space for studying nanofluidics. For a high- 
frequency nanomechanical resonator with resonance frequency u/2it, one can tune ujt over 
a wide range — in fact, possibly reaching the limits of the Newtonian approximation in a 
given liquid or gas. This not only allows experimental probing of a flow regime that was 
inaccessible by past experiments [71, |8( , but also presents the unique prospect of designing 
nanodevices for key technological applications. 
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To complement the recent theoretical interest in high frequency nanofluidics [2l, [9|, [10|, 111] , 



we experimentally studied the interaction of high-frequency nanomechanical resonators with 
a gaseous environment. The gaseous environment presents an ideal fluid for these studies, 
where one can effectively tune r by changing the pressure p. On the other hand, varying the 
resonator dimensions changes the mechanical resonance frequency uj/2it. When combined, 
the two experimental parameters allow ujt to be varied over several orders of magnitude 
effectively. 

In order to cover a broad frequency range, we fabricated silicon doubly-clamped beam 
resonators with varying dimensions w x hx I displayed in Table [I] using standard techniques 



12j. To further extend the frequency range, we employed fundamental and first harmonic 



modes of two commercial silicon AFM cantilevers (Tabled]). For the measurements, we used 
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TABLE I: Device parameters, transition pressure p(Wi = 1) and the approximate lower pressure 
limit Pmin for accurate measurements for the devices used in the study. 1 st harmonic mode was 
also employed for some AFM cantilevers. 
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a pressure-controlled optical characterization chamber connected to a high purity N2 source. 
We actuated the out-of-plane modes of the resonators electrostatically and measured the 
displacements optically [l2j. All measurements were performed under linear drive; more- 
over, the results remained independent of the rms displacement amplitudes of 1-10 nm as 
confirmed by Michelson interferometry. Figure [T] depicts the typical resonant response of a 
nanomechanical resonator as the background N 2 pressure in the chamber is increased. The 
frequency shift is due to the mass loading from the boundary layer j^, while the broaden- 
ing results from the energy dissipation in the fluid. The analysis can be simplified by using 
a one-dimensional damped harmonic oscillator approximation [ljj, x + ^x + uj 2 x = f /m, 
where f/m represents the force per unit effective mass of the resonator. The quality fac- 
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tor Q, which is a comparison of the stored energy to the dissipated energy per cycle, is 
related to 7 as 7 ~ lu/Q. Here, we extracted both the resonance frequency u/2n and Q 
using nonlinear least squares fits to the Lorentzian response of the resonator. In addition, 
for low-Q (high pressure), we verified the Lorentzian fit results through fits to the real and 
imaginary components of the complex transmission [14J]. Typical changes in u and Q of a 
nanomechanical resonator during a pressure sweep are shown in the inset of Fig. [TJ Both uj 
and Q approach their respective intrinsic values, ujq and Qo, at low pressure. 

Before presenting further results, we must clarify the nature of the fluidic energy dissipa- 
tion. The motion of the fluid with respect to the solid boundary creates a complex, position- 
dependent shear stress on the resonator surface. The inertial and dissipative components 
of the net shear force are proportional to the displacement and the velocity, respectively. 
For a single device, as the pressure is changed, the fluidic dissipation can be quantified by 
either the fluidic quality factor Qf given by Qj 1 = Q^ 1 — Qq 1 or the fluidic dissipation 
jf = u/Qf. To compare different devices with varying sizes and geometries, one needs to 
further realize that the fluidic dissipation is proportional to the effective surface area S e ff, 
while the stored energy in the resonator is proportional to the effective mass m 24j]. With 
this naive assumption, we define a normalized fluidic dissipation, 7 n = / ~ffm/S e ff. The lower 
pressure limit p m in for accurate Qf measurement is set by Qo'- as one approaches p m i n , the 
intrinsic losses in the resonator dominate the measurement. The upper limit is 1000 Torr. 
Table fl] displays p m in for each device along with Q . 

The normalized fluidic dissipation 7„ observed in three different resonators are presented 
in Fig. El^a),(b) and (c) as a function of gas pressure. A change in the slope is noticeable 
for the data in (a), and (b) at approximate pressures of 1 and 300 Torr, respectively. The 
turn points marked in the plots correspond to lot w 1, and is discussed in detail below. 
For the highest frequency beam at 102.5 MHz shown in (c), the turn point falls outside 
the available pressure range, i.e., 1000 Torr. Molecular flow model 10], which takes into 
account specular collisions, fits our data only at the ideal gas limit at low pressure. Note 
ihat a multiplicative constant of 0.9 was used in all three to improve the fits. Viscous effects 



m 



16| and squeeze-film effects 17|, ll8[], commonly observed in MEMS, did not introduce 



significant damping for the small high-frequency devices up to atmospheric pressure. 

The solid lines in Fig. [2] are fits to a theory by Yakhot and Colosqui [9j developed from the 
Boltzmann Equation in the relaxation time approximation. After geometric normalization 
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FIG. 1: Resonance of a silicon doubly-clamped beam of width w = 500 nm, thickness h = 280 
nm and length I = 11.2 fim at various N2 pressures in the chamber: p = 0.049, 5.4, 32, 100, 302 
and 942 Torr. The inset shows the extracted quality factor Q and normalized resonance frequency 
uj/ujq of the same device as a function of pressure. Here, Qq « 1530. 




Pressure (Torr) 



FIG. 2: Normalized fluidic dissipation y n as a function of pressure for (a) a cantilever with dimen- 
sions (id x /i x I) 53 x 2 x 460 /im, and nanomechanical doubly-clamped beams with dimensions of 
(b) 230 nmx200 nmx9.6 /mi (c) 240 nmx200 nmx3.6 /im. Resonance frequencies are as indicated 
and the approximate turning points are marked with arrows. The lines are fits to molecular colli- 
sion model [l(| and Eq. [T] using r ~ 1850/p [9J. The molecular collision and the Yakhot-Colosqui 
predictions were multiplied by 0.9 and 2.8, respectively, for all resonators, (d) Relaxation time r as 
a function of pressure. The points were extracted from fluidic dissipation data sets, such as those 
shown in (a) and (b), of 13 resonators. The solid line is a least-mean-squares fit and indicates that 
r w 1850/p. 
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FIG. 3: (a) Normalized fluidic dissipation ~f n as a function of the resonator frequency loq for several 
resonators at four different pressures. From top to bottom, r ~2.3, 4.6, 9.2, 18.5 ns. The lines 
are fits calculated using Eq. [TJ Wi ~ 1 points are marked with an arrow for each pressure. 
(b)The scaling of j n / ^Juippf for all resonators with lot. Each symbol corresponds to an individual 
resonator. All the predictions were multiplied by the same fitting factor of 2.8 (also see Fig. [2]). 

and imposing the no-slip boundary condition, this theory culminates in the expression 
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COS 



'l+UJ 2 T 2 ) 3/i 



tan ujt 



[(1 + wt) 



/ tan lut \ 

In Eq. [TJ 7 n is expressed in terms of the viscosity /i, the density pj, the effective relaxation 
time r of the fluid, and the frequency uj of the resonator. The only unavailable parameter is 
the relaxation time r. In order to obtain the fits, we assumed that r satisfied the empirical 
form, r cx 1/p 0, Q. The key prediction of the Yakhot-Colosqui [9] theory is that the turn 
point in 7„ occurs when r ~ l/u. Thus our experiments provided a direct and unique way 
to extract r as a function of pressure p: Fig. EJ^d) displays the experimentally extracted r 
from the transition points of multiple resonators as a function of pressure. Through linear 



fitting, one can obtain the expression r »s 1850/p [in units of nanoseconds when p is in 
units of Torr]. The end result of this exercise are the self-consistent fits in Fig. |2](a), (b) 
and (c). To improve the fits, the results emerging from Eq. [JJwith the appropriate material 
properties and r were multiplied by 2.8. In general, all our data sets could be fit adequately 
using Eq. [TJ after multiplying by 2.8 ± 0.7. 

The fluidic dissipation in individual resonators shown in Fig. [2] suggests that there, in- 
deed, is a transition at lot « 1, obtained by tuning r. Further support for this transition 
comes from extended measurements in the frequency parameter space. Figure EJ^a) shows 
7 n from different resonators spanning a huge frequency range. Here, 7„ is plotted against 
the resonator frequency ujq at four different pressures, i.e., four different r. This comparison 
between different devices with different sizes is possible only after normalization of the dis- 
sipation by S e ff/m [25]. The solid lines in Fig. [3](a) are fits to Eq. [fusing t ~ 1850/p. The 
points marked by arrows correspond to lot ~ 1. Again, we have multiplied all the fits by 2.8 
as in Fig. [2l This multiplicative constant probably arises from adapting the theoretical ex- 
pressions 9[ for an infinite plate oscillating in-plane to the finite and rectangular resonators 
oscillating out-of-plane. Our surface-to-volume normalization does not give the absolute 
dissipation, while it appears to be useful for comparing different devices. The elucidation of 
finite size effects in complex geometries is the subject of our ongoing computational research. 
Figure [3](b) shows all the data, 7„/ y/ujjlpj, from all devices collapsed onto a single curve, 
plotted against over four decades of the dimensionless parameter Wi = lot. Each symbol 
type in Fig. [3](b) corresponds to a separate resonator; the Wi is obtained by multiplying the 
resonator frequency by the corresponding r from r « 1850/p. The solid curve is obtained 
from Eq. [TJ 

Also apparent in Fig. [TJ is a small decrease in the resonance frequency as the pressure is 
increased. The observed decrease is primarily due to the mass of the fluid mf that is being 
displaced in-phase with the resonator: ~ ^ ^ [4]. Both the geometry and the frequency 
of the moving surface is expected to play a role in determining m./ and thus Aoo. In an effort 
to rule out the geometry effects, we studied the frequency shift in four beams of identical 
widths (w = 500 nm) and thicknesses (h = 280 nm) but varying lengths and resonance 
frequencies as a function of pressure (Fig. @|. Here, the plotted Alo/lo corresponds to the 
approximate fluid mass per unit beam length. The arrows mark Wi ~ 1 for each beam and 
the dashed lines represent best line fits corresponding to Alo oc p 1 ^. The molecular flow 
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FIG. 4: Normalized frequency ^ for beams of identical widths w = 500 nm and thicknesses 
h = 280 nm, but varying lengths and resonance frequencies. The pressure for Wi ~ 1 is marked 
with an arrow . Dashed lines are least-mean-squares fits to Aoj oc p 1//3 . 

model [loj], which is appropriate at low pressure, does not predict any mass loading and 
frequency shift. At high pressure, the Stokes' expression for an oscillating plate [l[ can be 
used to obtain an approximation for the boundary layer thickness, 5 = ^jjzi- This, however, 
results in a pressure independent Aa>, given that typical 5 ~ 1 — 5 /im, and consequently, 



~m ~ p^/r Yakhot-Colosqui theory [9| underestimates the magnitude of the observed 
frequency shift. We expect to understand the nature of the boundary layer in the near future 
by studying the scaling of experimental frequency shifts in a wide range of geometries and 
frequencies, and through supporting computational analyses. 

The transition observed in our experiments can be interpreted in most general terms 
as follows. The simple linear relation between stress and rate-of-strain in a Newtonian 
fluid breaks down at high frequencies. The Boltzmannian theory developed by Yakhot and 
Colosqui suggests that this result is independent of the nature of the fluid. This transition at 
Wi = lot « 1 was described [9j as a "viscoelastic to elastic" transition owing to the fact that 
the waves generated in the fluid by the resonator motion do not decay as uj —>■ oo. There 
also appears to be some universality with respect to device geometry: both in cantilevers 
and doubly-clamped beams, the same naive geometric normalization resulted in a consistent 
analysis. 

There is a relentless effort to develop nanomechanical resonators operating in gaseous 
[20I and liquid environments 21]. Our results should impact the design of next-generation 
nanomechanical resonators. Figure [3] suggests that fluidic dissipation saturates at high 
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frequencies. Take, for instance, two doubly-clamped beam resonators with identical widths 
and thicknesses but different lengths, i.e., identical ~ — + \ but different frequencies 
such that <jj\ < 10%. If U\ < 002 < 1/t, the ratio of the quality factors of the two resonators 
in fluid becomes ~ . — . On the other hand, if oh < 1/r < cu 2 , then S^i ^ w 
Finally, for 1/r < U\ < ~ Thus, shorter, higher-frequency resonator will always 

be more resilient in a given fluid but the degree of resilience depends upon the fluid r. Yet, 
for two devices with identical frequencies, the smaller one with the larger will have the 
lower Qf. For the case where both and device frequency increase, the nature of the 
scaling determines the end result. Finally, the surface roughness, especially for very small 
devices, is expected to have an important role in nanofluidics of nanomechanical resonators 



22]. 
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